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We study the effects of metallic doping on the electronic properties of graphene using density 
functional theory in the local density approximation in the presence of a local charging energy 
(LDA+U). The electronic properties are sensitive to whether graphene is doped with alkali or 
transition metals. We estimate the the charge transfer from a single layer of Potassium on top of 
graphene in terms of the local charging energy of the graphene sheet. The coating of graphene with 
a non-magnettc layer of Palladium, on the other hand, can lead to a magnetic instability in coated 
graphene due to the hybridization between the transition-metal and the carbon orbitals. 
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PACS numbers: 71.15.Mb, 73.20.At, 75.10.Lp 



INTRODUCTION 



Graphene, a two dimensional (2D) allotrope of car- 
bon on a honeycomb lattice, is characterized by el- 
ementary electronic excitations described in terms of 
Dirac fermions^, a solid state realization of a relativis- 
tic system^. The observation of theoretically predicted^"^ 
anomalous plateaus in the quantum Hall effecli^i^ and of 
a universal minimum of conductivity^ attracted a lot of 
interest. In the absence of doping, graphene behaves as 
an unusual metal with low density of states^^. Because 
the linear band spectrum is a robust feature of the hon- 
eycomb lattice, the excitations in graphene are particles 
with zero effective mass that propagate coherently very 
large distances disregarding the amount of impurities, al- 
lowing graphene to sustain supercurrents- . By applying 
a bias voltage, the carrier density of graphene can be con- 
trolled by electric field effect allowing for many practical 
applications ranging from the production of electronic 
lenses^ to the fabrication of semiconductors with a tun- 
able gap in bilayeraiS. 

One of the roads still unexplored in the material sci- 
ence of graphene is the tailoring of its electronic prop- 
erties by chemically adsorption of metallic atoms. The 
electronic properties that result from adsorption depend 
strongly on the ionic and/or covalent character of the 
bonds formed between carbon and the metal. Alkaline 
metals are good electron donors because of the strong 
ionic character of their bonding, increasing dramatically 
the number of charge carriers in graphene^. At the 
same time, the metallic bands can be strongly affected 
by the presence of the graphene lattice, as observed in 
some graphite intercalated compounds (GIC)^^. On the- 
oretical grounds, both effects can produce superconduc- 
tivity in coated graphene^^. Transition metals by their 
turn tend to make strong covalent bonds. Because of the 
strong electron-electron interactions in these materials, 
they are more susceptible to induce magnetic instabili- 
ties. The tuning of magnetism in coated graphene can 
open other routes to spintronics through new spin- valve 
devices'^ . Similar effects have been studied in the con- 
text of carbon nanotubes^-. 



We study the effect of metallic coating of graphene 
using two different metal atoms, potassium (K), and pal- 
ladium (Pd). Among alkaline metals, K atoms have a 
particularly good matching with the graphene lattice and 
their adsorption on graphite surfaces has been a topic of 
extensive researcbi^. We address the problem of charge 
transfer between K and graphene, and discuss the nature 
of the K conduction band in coated graphene. Although 
there are many transition atoms that interact strongly 
with carbonic, Pd has very polarizable bands and hence, 
is a natural candidate for the generation of magnetism in 
coated graphene. 



II. NUMERICAL CALCULATION 

We investigate the electronic properties of coated 
graphene from a band structure calculation based in den- 
sity functional theory (DFT) in local density approxima- 
tion (LDA). In order to study a single-layer of graphene, 
we repeat graphene layers periodically with 33.6 A of 
vacuum separation. The electronic structure is calcu- 
lated using the all-electron full-potential linearized aug- 
mented plane wave (FLAPW) method" with corrections 
to the exchange-correlation potential calculated in the 
generalized gradient approximation (GGA)^*. To fur- 
ther take into account the local interactions, we go be- 
yond GGA assigning a phenomenological parameter U 
(LDA-|-U) ^^i^° , which corresponds to the effective local 
potential between two electrons in the same localized or- 
bital. 



A. K coated graphene 

In K coated graphene, K atoms sit on top of the 
hexagons around 2.7 A away from the graphene planc^^.. 
Since K is much larger than carbon, the maximum cov- 
erage is achieved at one K atom per 8 carbons, KCs, 
what corresponds to one monolayer. In this concentra- 
tion, we relax the K-graphene distance to 2.68 A within 
GGA after the minimization of the interatomic forces is 
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Figure 1: (color online) On the left: superposition of the 
bands for one monolayer of K (red and brown lines) and pure 
graphene (black lines) in the unit cell of KCg. On the right; 
band plot of KCg in GGA ((7 = 0). 



done. Since the C-C bonds are much stronger than the 
bonds of C with the metal, the effect of the relaxation 
of the graphene lattice due to the coating is a very small 
effect^^. We use the lattice parameter of pure graphene, 
a — 1.42 A, in the calculations. 

When K is deposited on top of graphene, part of the 
electrons in the 4s-band are transfered to the 7r*-band of 
C to compensate for the difference in electronegativities. 
In this process, the 7r-bands suffer a rigid shift, generat- 
ing a pocket of electrons with a finite density of charge 
carriers in the graphene plane. At the same time, there is 
a significant redistribution of charges within the K states 
which accompanies the formation of K-C bond o^^'^'^ . In 
the plot on Fig. 1 (left), we show the bands of graphene 
in the expanded unit cell of KCg superimposed with the 
bands of a free standing K monolayer. The 4s band of K 
shown in red forms a nearly circular Fermi surface around 
the center of the Brillouin zone (BZ) at the F point, while 
the TT-bands of pure graphene cross the Fermi level at 
the corners of the BZ (K point). After a careful com- 
parison between the high energy bands of KCs with the 
K and graphene bands in separate, we are able to find 
a nearly one-to-one correspondence in the energy range 
0-6 eV between the empty bands of KCs, shown in Fig. 
1 (right), and the empty bands of the free K monolayer. 
However, we find no trace of the graphene nearly free 
electron (NFE) bands in KCg, which in pure graphene 
start from ~ 3 eV above the Fermi level, as shown in 
Fig. 1—. This suggests that the interstitial NFE bands 
drop to the Fermi level and hybridize with the 4s band 
of K in a similar way to previous band structure analysis 
on GIC^^ and carbon nanotubes^-. The strong downshift 
of these high energy bands is an electrostatic response of 
the graphene NFE bands to the potential induced by the 
background of positive charge of the K ions^S. This pic- 
ture indicates that the charge of this system is distributed 
in the graphene layer and also in the interstitial space of 



the KCs bilayer. 

The amount of charge transfer to graphene in KCg is 
still an open question in the literature and has been stud- 
ied for many years in K deposited graphite surfaces"'^^. 
Previous LDA calculations with specific assumptions on 
the pseudopotentials have predicted a charge transfer of 
^ 0.17 e per K^^'^^ on KCs and ~ 0.1 e per K for one 
monolayer of K on top o a stack of graphite layers^LS^. 
On the other hand, a recent all-electron GGA cluster 
method calculation on KCs has predicted a substantially 
larger charge transfer of 0.46 e per K2S. We remark how- 
ever that because C orbitals are poorly screened by K, 
the local Hubbard potential in the carbon pz orbitals is 
an important input which modifies the difference of elec- 
tronegativities between graphene and the metal. Indeed, 
we observe within the LDA-I-U framework that the total 
amount of charge transfer is very sensitive to U. To esti- 
mate the charge transfer, we extract the energy shift of 
the K 4s band (red curve in Fig. 1) when the K layer 
is deposited on top of graphene. For U = 0, our DFT 
calculation predicts a charge transfer of 5Q = 0.51 e per 
K. As U increases the charge transfer is linearly sup- 
pressed and eventually extrapolates to zero for U ^ 25 
eV (see Fig. 2). The hopping energy of the electrons in 
graphene is i ~ 3eV. As in fuUerenes^S, if a Hubbard 
U ^ it can be taken as a reasonable estimate for the 
order of magnitude of U in graphene, then the LDA-I-U 
calculation may lower considerably the amount of charge 
transfer predicted by a pure GGA calculation^^. On the 
other hand, since the LDA-I-U method does not take po- 
larization effects into account, the local charging energy 
U can be additionally screened by the pocket of electrons 
in graphene. Our effective ?7 is a screened local charging 
energy in the graphene layer. 

The electronic density transfered to the graphene layer 
is 5Q/Ay^ where — 6\/3a^ is the area of the K unit 
cell on top of graphene. For a small charge transfer 5Q, 
the shift in the chemical potential of the 7r-bands is 

5iM= {vFla)\jTT/{<oV?,)5Q -2.3V^eV (1) 
where vp = 6eV A is the 7r*-band velocity around the K 
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Figure 2: Charge transfer SQ from K to graphene in KCg (in 
electrons per K) as a function of the effective Hubbard U in 
the C Pz orbitals. The charge transfer for [/ = has been 
calculated in GGA, for comparison with Ref.— . 
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point. Hence, we see that a very small charge transfer of 
6Q ~ 0.01 e per K is already enough to exceed the max- 
imum chemical shift achievable by the application of a 
bias voltage in graphene (« 0.2 eV^). For a charge trans- 
fer larger than ^ 0.2 e per K, the shift of the vr-bands is 
above 1 eV and reaches 6fi ^ 1.2 eV for U = 0, when 
the charge transfer is maximum (see Fig. 1). A 7r-band 
shift of the same order has been very recently measured 
with angle resolved spectroscopy (ARPES) for K coated 
graphene (KCg) on a SiC substrate^. In this experi- 
ment, the charge transfered from the substrate shifts the 
TT bands of pure graphene in 0.45 eV, which corresponds 
to ~ 0.04 e per unit cell of KCs- The total 7r-band shift 
was estimated in ~ 1.2 eV, indicating that each K atom 
transfers the order of 0.5 e to graphene. This result con- 
trasts with a recent photoemission experiment for a K 
monolayer adsorbed on the surface of graphite, which es- 
timated the charge transfer in ~ 0.1 e per K'^^. According 
to our DFT results, these experiments suggest that the 
local charging energy of graphene is strongly screened by 
the presence of the SiC substrate (see Fig. 2). 

In K intercalated graphite, previous DFT pseudopo- 
tential calculations have reported charge transfers larger 
than 0.7 e per 

K33^^ These results may indicate 
that the character of the K/C bonding is modified by 
intercalation'^'*. A more conclusive comparison between 
GIC with the present results, however, requires a closer 
examination of the intercalated case through a similar 
DFT method. 

The size control of the charge pockets in graphene can 
be achieved by dilution of the K coverage. For dilution 
less than 0.9 of a monolayer, however, K does not form 
a uniform metallic lattice due to clustering, leading to 
insulating behavior—. 
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B. Pd coating case 

Because of the low coordination number, the d orbitals 
of transition metals can be very localized in systems of 
low dimensionality. Among 4c? transition metals, bulk 
paramagnetic elements like Ru, Rh, and Pd exhibit mag- 
netism on surfaces and in nanosystem a'^^i'^^'^^ . In partic- 
ular, Pd atoms are not intrinsically magnetic because of 
the closed shell configuration of the 4c? orbitals (4d*°5s°). 
In bulk, the s and c? bands of Pd hybridize (4c?*°~'^5s'^), 
and Pd exhibits strong Pauli paramagnetism with a high 
magnetic susceptibility. Due to the high density of states 
near the Fermi surface, bulk fee Pd is close to a ferromag- 
netic instability. Expanding the fee lattice in 5%, the 
Stoner criterion is satisfied and Pd becomes an itinerant 
ferromagnet due to the enhancement of the s — d band 
hybridization?^^. 

On top of graphite, Pd atoms do not form a uniform 
metallic coating but large clusters^. It is not clear to 
what extent Pd atoms can grow laterally to form islands 
as Ru^S, and Rh"*" or if they agglomerate to form clusters 
of a few layers thick. If we assume that the Pd atoms 



Figure 3: (color online) Band plots for (a) graphene, (b) a 
free standing monolayer of Pd for [7 = and (c) free Pd 
monolayer for U = 11 eV. In plots (d), (e) and (f) we display 
the bands of PdC2 for: (d) U = 0, (e) U =11 eV in Pd for 
spin down and (f) U = lleV in Pd for spin up. 



inside the islands form a low temperature commensurate 
structure with the graphene lattice with the same period- 
icity as in Rh adsorbed on the surface of graphite^°, the 
Pd atoms sit at the center of the hexagons of graphene 
separated by a distance of V3a ~ 2.46A from the next Pd 
neighbor. In this case, the Pd-graphene equilibrium dis- 
tance is 3.35 A, after we minimize the interatomic forces 
between the Pd and graphene layers. This configuration 
represents a lateral compression of ^--^7% with respect to 
the bond length of Pd in a square lattice^, and a very 
small magnetization is expected as a result of the weak 
hybridization of the bands^. 

Nevertheless, the splitting of the bands due to the hy- 
bridization of Pd and graphene orbitals at finite U can 
produce a strong enhancement of the density of states at 
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Figure 4: Magnetization of Pd coated graphene, PdC2, (in 
Bohr magneton) vs. local potential U in the Pd orbitals. 



charge neutrality and the redistribution of charge among 
the hybridized d orbitals due to the local charging energy 
enhances dramatically the density of states in the local- 
ized bands of Pd at the Fermi surface. In this situation, 
the Pd bands strongly polarize, as displayed in Fig. 3 
(e) and (f), and we observe a total magnetization of 0.64 
/is in the Pd unit cell (see Fig. 4), indicating a ferro- 
magnetic instability. The strength of the magnetization 
is very sensitive to the specific value of the assigned local 
potential for t/ > 6 eV, where the Stoner condition is 
satisfied and the system has a ferromagnetic phase tran- 
sition, as shown in Fig. 4. 



the Fermi surface. We show that the Pd bands in coated 
graphene can strongly polarize through a Stoner transi- 
tion. We focus in the Pd coating to illustrate the effect of 
the metal-carbon hybridization in the production of itin- 
erant ferromagnetism. Since the orbitals in the transition 
metal monolayer are more localized than in graphene, we 
assign in first approximation a single effective parameter 
U to the metallic bands in our LDA+U calculation. 

In Fig. 3(b) we show the spectrum of a free standing 
monolayer of Pd in the unit cell of graphene. The bands 
in red can be associated at the F point (the center of the 
BZ) to the two degenerated bands and 4c?j,2 of the 
free Pd atom. The two bands that cross the F point at 
^ —3 eV can be associated near F to the degenerated 
bands 4:dxy and 4d^2_j^2 of free Pd, while the 5s-band 
crosses the F point at ^ — 2eV. Finally, the high energy 
band which crosses the F point at ~ 6 eV below the Fermi 
level corresponds to the 4(i^2 orbital of free Pd around the 
center of the BZ. 

When we turn on the local potential to U = 11 eV 
in the free Pd monolayer [see Fig. 3(c)], we clearly see 
that the bands in red suffer a rigid downshift of ^ 1 
eV, indicating the presence of localized states in these 
bands. The band which is nearly ffat between points K 
and M in Fig. 3(b) shifts in ^ 1 eV around the bor- 
der of the BZ (points K and M), where the effects of the 
lattice symmetry are stronger, and remains unaltered at 
the center of the BZ, where the d^-2_y2, dxy orbital sym- 
metries dominate, suggesting that the states around F in 
this band are more delocalized along the plane of the Pd 
monolayer. The other two Ad bands and the 5s band are 
almost insensitive to the local potential and are delocal- 
ized bands. We do not see any trace of magnetization 
in the free standing Pd monolayer for U in the range of 
0-12 eV. 

In the sequence of plots, in Fig. 3(d), (e) and (f), we 
put the Pd monolayer on top of graphene for [/ = and 
then we set U finite. Because graphene is more elec- 
tronegative than Pd, the tt band is shifted up in ~ 1 eV 
for U = Q. In this case, the d orbitals of Pd strongly 
hybridize with the graphene Pz orbitals but the Stoner 
condition is not satisfied and the system does not magne- 
tize. For [7 = 11 eV, however, the 7r-bands nearly recover 



III. SUMMARY AND CONCLUSION 

We have explored the electronic properties of graphene 
coated with alkaline and transition metals from a 
band structure calculation in LDA+U. We propose that 
graphene can be tailored either by the control of the num- 
ber of charge carriers, which can be tuned by coating 
graphene with a controlled coverage of alkaline metals, 
such as K, or by inducing a ferromagnetic instability in 
graphene through the coating with a transition metal. 

We addressed the problem of the charge transfer be- 
tween one monolayer of K and graphene, emphasizing 
the role of the local charging energy in the graphene 
layer. Our DFT results in coated graphene suggest that 
the metallic band of K hybridizes with the carbon NFE 
band in a similar way as previously reported in GIG and 
carbon-nanotubes. We emphasize that the problem of 
coating is conceptually different from the adsorption of 
isolated alkaline adatoms in graphene, where the metallic 
orbitals strongly localize around the impurities. In par- 
ticular, the estimation of the charge transfer from each 
K impurity in the dilute phase depends on specific defi- 
nitions of the atomic charge^. In graphene coated with 
a uniform monolayer of K, the K s electrons delocalize 
and the charge transfer can be accurately estimated by 
the energy shift of the bands at the Fermi level, as shown 
in Fig. 1. 

The amount of charge transfered from K within the 
LDA-t-U framework is sensitive to the effective local 
charging energy assigned to the carbon orbitals. For 
t/ = 0, we find a maximum charge transfer of 0.51 e per 
K, which is linearly suppressed with the increase of the 
effective charging energy U . A comparison of our DFT 
calculation with recent spectroscopy and photoemission 
experiments suggest that different substrates can give rise 
to very different values of U in graphene. 

For Ad transition metals such as Pd, we showed that 
the hybridization of the graphene Pz orbital with the lo- 
calized d orbitals at finite U can produce strong itiner- 
ant magnetism in coated graphene. The onset of ferro- 
magnetism in Pd coated graphene is identified through a 
Stoner transition for ?7 > 6 eV, where the local bands of 
Pd polarize. 
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